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Abstract: This is the second of a series of three papers examining how viable it is for 
entanglement to be sustained at high temperatures for quantum systems in thermal equilib¬ 
rium (Case A), in nonequilibrium (Case B) and in nonequilibrium steady state conditions 
(Case C). The system we analyze here consists of two coupled quantum harmonic oscil¬ 
lators each interacting with its own bath described by a scalar held, set at temperatures 
Ti > T 2 . For constant bilinear inter-os dilator coupling studied here (Case Cl) owing to 
the Gaussian nature, the problem can be solved exactly at arbitrary temperatures even for 
strong coupling. We hnd that the valid entanglement criterion in general is not a function 
of the bath temperature difference, in contrast to thermal transport in the same NESS 
setting [1]. Thus lowering the temperature of one of the thermal baths does not necessarily 
help to safeguard the entanglement between the oscillators. Indeed, quantum entangle¬ 
ment will disappear if any one of the thermal baths has a temperature higher than the 
critical temperature Tc- With the Langevin equations derived we give a full display of how 
entanglement dynamics in this system depends on Ti, T 2 , the inter-oscillator coupling and 
the system-bath coupling strengths. For weak oscillator-bath coupling the critical temper¬ 
ature Tc is about the order of the inverse oscillator frequency, but for strong oscillator-bath 
coupling it will depend on the bath cutoff frequency. We conclude that in most realistic 
circumstances, for bosonic systems in NESS with constant bilinear coupling, ‘hot entan¬ 
glement’ is largely a fiction. In Paper III we will examine the case (C2) of time-dependent 
driven coupling which contains the parametric pumping type described in [2] wherein 
entanglement was first shown to sustain at high temperatures. 
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1 Introduction 

Recently Galve et al [2] (see also [3]) pointed out the possibility of keeping quantum en¬ 
tanglement alive in a system at high temperatures by driving the system of two oscillators 
with a time-dependent interaction term. This has generated a great deal of interest in 
understanding the underlying issues and the basic mechanisms of obtaining the so-called 
‘hot entanglement’ [4]). The word ‘hot’ conveys three layers of meaning in three different 
contexts, referring to quantum systems A) kept in thermal equilibrium at all times, B) in 
a nonequilibrium condition and evolving and C) in a nonequilibrium steady state at late 
times. Thus before making sweeping statements one needs to discern and analyze systems 
under at least these three separate situations for the behavior of quantum entanglement 
therein. 

We have analyzed Case B) described above in some detail in our first paper, obtaining 
the parameter ranges for entanglement to survive at a hnite temperature and comparing 
with the results for Case A obtained earlier in e.g., [5, 6 ]. Our results indicate that, gener- 
ically, when two coupled oscillators separated at a fixed distance evolve under the influence 
of a shared thermal bath, their dynamics is usually highly non-Markovian. The asymptotic 
correlation /entanglement between the oscillators tends to survive better under 1 ) stronger 
inter-oscillator coupling, 2) weaker oscillator-bath interaction and at 3) a shorter distance 
between them. In the case of weak oscillator-bath coupling, the critical temperature is still 
bounded by the inverse oscillator’s natural frequency, but tends to be lower than that the 
critical temperature in Case C, due to the hnite separation between the oscillators. The 
largest separation before the entanglement drops signihcantly is of the order of the inverse 
cutoff frequency inherited in the thermal bath, and the distance will decrease with higher 
bath temperature. In this limit, the results are similar to Case A. This is not unexpected 
since it is known [26, 27] that in the weak coupling limit, both conhgurations will yield 
similar results; furthermore, the non-Markovian mutual interaction between the oscillators 
is minimal in the weak oscillator-bath coupling regime. For stronger oscillator-bath inter¬ 
action, the mutual interaction can sustain over a very long history in the evolution of both 
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oscillators. Deviation in results between Case A and Case B will emerge. Nonetheless a 
strong oscillator-bath interaction can likely induce dynamical instability in the oscillators, 
a case worthy of closer analysis later. 

In this paper we analyze condition C) where the system can maintain a nonequilib¬ 
rium steady state (NESS) at late times. Since NESS is a distinctly generic state, playing 
an important role for nonequilibrium systems as fundamental as the equilibrium state in 
quantum statistical mechanics, it is important to clarify the behavior of high temperature 
quantum entanglement under such conditions. We illustrate these two conditions with two 
generic models: Case B) is exemplified by a quantum system made of at least two har¬ 
monic oscillators (HO) interacting with one common thermal bath] Case C) is exemplified 
by a quantum system composed of two coupled harmonic oscillators each interacting with 
its own (private) thermal hath. We wish to inquire about how entanglement initially present 
between the two quantum oscillators evolves in time, and calculate at what temperature 
(approaching from below) it begins to die out. 

To identify the root cause of quantum entanglement existing at high temperatures, 
if it does at all, one needs to identify the determining factors. Coupling in the system 
is certainly an important factor. Intuitively the stronger the coupling in the system, the 
weaker the coupling of the system to the baths, the better preserved the entanglement will 
be. If the coupling can be tuned to “cruise alongside” how entanglement evolves in time, 
to even amplify it along the way, the better the chance of keeping the entanglement alive. 
To see these effects more clearly we further divide the nonequilibrium steady state cases 
into two subcases. Cl and C2. Case Cl is for time-independent inter-oscillator coupling, 
and Case C2 for time-dependent inter-oscillator coupling. Before one can bring these 
cases under the same roof of nonequilibrium steady state condiiton one needs to prove or 
demonstrate that indeed a steady state exists at late times in these setups. We have so far 
shown the existence of NESS only for Case Cl in [1]. ^ . 

Before we treat the Case Cl scenario in full which is the main goal of this paper, we 
first give a brief description of a Case C2 model to mark the differences so the results of 
our work can be placed in perspective. 

As a model for Case C2 the system is made up of two quantum oscillators interacting 
with each other via a time-dependent (sinusoidal) coupling. Unlike Case Cl where the 
temperatures of the two baths are different, here they could be the same. In fact the 
temperature of the thermal bath and how strong the oscillators are coupled to the baths 
are not important. The nonequilibrium condition is provided by the external driving agent. 
Driving leads to production of entanglement even at very high temperatures. Eor instance, 
even with a weak environmental coupling, a strong driving amplitude still provides a higher 
critical temperature. 

The physics for these two cases albeit both in NESS is also very different. As explained 

^It naturally behooves upon advocates of hot entanglement [4] under NESS, namely, those with time- 
dependent coupling as exemplified by [2] to prove or show the existence of a NESS under those conditions. 
It may not be a straightforward task. In fact, for lack of a proof that systems with time-dependent coupling 
can approach NESS it is probably more prudent to call this setup Class D, and only after such a proof shall 
one reinstate it into Class C for systems which admit NESS. 
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in [2] , it is the squeezing of the system provided by the external agent and the parametric 
amplihcation (pumping) which can offset the thermalization /equilibration process natu¬ 
rally expected for the systems interacting with a bath and dominant at high temperatures. 
Parametric driving is what sustains the entanglement in the system. We will study this 
case in our sequel paper. 

1.1 Time-independent bilinear inter-oscillator coupling 

In the case of a chain of quantum harmonic oscillators coupled bilinearly with each other 
and with the baths the dynamics of the total system admits a complete solution, by virtue 
of its Gaussian nature, for all temperatures and for strong coupling within the system and 
with the baths. This model has been studied by many authors [10-12, 14, 15]. In our recent 
work [1] functional methods are used to provide an explicit demonstration of the existence 
of a nonequilibrium steady state. Here we apply the results obtained therein to a study 
of quantum entanglement in NESS, with the aim of quantifying the claims made in the 
literature [4], alluding to the possibility of entanglement survival at high temperatures for 
systems in NESS. Note the present setup of bilinear coupling is different from that of Galve 
et al [2] where the interaction between the two oscillators is via parametric pumping. Eor 
this setup a recent paper closest to our intent is that by Ghesquire, Sinayskiy &: Petruccione 

[ 7 ]. 

1.2 Comments on Claims by Other Authors 

We make a brief summary of what GSP have done and what claims they made below. 

Eor the same model as mentioned above, namely, two bilinearly coupled quantum 
harmonic oscillators each interacting with its own bath GSP derived a perturbative ‘pre- 
Lindblad’ master equation without invoking the rotating wave approximation (non-rotating- 
wave, or NRW) [8]. They consider two situations: Eor the study of entanglement they 
consider the high temperature regime in their Eq. (3) valid for both strong and weak 
interaction strength with the baths. Eor the consideration of entropy dynamics related to 
equilibration issues they take the weak system-bath coupling limit and arrived at their Eq. 
(4). We will only be concerned with the entanglement issue here. GSP made the following 
claims: 

a) Entanglement persists for longer times at lower temperatures. 

b) In the weak system-bath coupling limit, the late time steady state developed is 
independent of the initial conditions. 

c) Eor the equilibrium case, there exists a critical temperature which is consistent with 
the result of [5] in the limit. 

We limit to two comments regarding their method and claims here. The major differ¬ 
ences will become clear in our results with quantitative representation via graphs found in 
later sections. 

1) Regarding the method and approximations: A perturbative ‘pre-Lindblad’ master 
equation, even without the RWA, does not in general satisfy the complete positivity condi¬ 
tion. Although it works better for strong coupling to the environment the results obtained 
under these approximations have unphysical behavior at low temperatures. Eor example. 
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Ludwig et al [12] pointed out the effect from the environment cutoff has to be handled with 
care. 

2) The claim statements are too general - they may not hold for specific conditions. 
They need be qualified more carefully by specifying the range of (in) validity of the ap¬ 
proximations introduced. E.g., Point a) above is sort of expected, but does it also imply 
that entanglement can be generated and be sustained if the temperature of both baths are 
sufficiently low, even though the system state is initially separable? Point b) regarding the 
existence of a NESS - it has been demonstrated for arbitrary strength in bilinear inter¬ 
oscillator coupling and for arbitrary temperatures of the two baths [1]. Point c) There is a 
distinction between i) a system of two coupled oscillators each with its private bath under 
NESS studied here, setting the two baths to be at the same temperature (presumably what 
their ’equilibrium’ condition entails) and ii) the system in one common thermal bath (what 
we call Case A). The situation is a lot more complex - see discussions in the last section 
of this paper. 

The above questions and a broader set of issues will be addressed in a fuller treatment 
of this generic (bilinear coupling) NESS model in the sections below. 

1.3 Our Method and Key Findings 

The model we use in this work to describe entanglement dynamics at high temperature, 
namely, two coupled oscillators each interacting with its private bath at different temper¬ 
atures, has been treated in full in our earlier paper [ 1 ], where one can find more technical 
details of the whole framework. Entanglement in a harmonic chain is also a well-explored 
subject . The Gaussian nature of this model allows us to obtain exact solutions for ar¬ 
bitrary coupling strengths and temperatures . The central quantity to calculate is the 
covariance matrix at finite temperature and at late times, where it has been shown that 
the system approaches NESS. The Peres-Horodecki-Simon entanglement criterion [16-18] 
can be calculated without approximation. This approach has been shown to be totally 
equivalent to that of directly deriving the reduced density operator of the system [ 1 , 19]. A 
short way to report on our hndings is that quantum entanglement will disappear when the 
bath temperatures become higher than a critical temperature [T^ = l//3c)- Also not sur¬ 
prisingly, asymptotic entanglement is easier to sustain for stronger inter-oscillator coupling 
and weaker oscillator-bath coupling. The true gain of this investigation is a full display 
via the Langevin equations we derived of the dependence of entanglement dynamics on 
the three parameters in this model, temperatures {Ti,T 2 ) of the baths, the intra-system 
(inter-oscillator) coupling a and the system-bath coupling strengths 7 . Their interplay is 
presented in the plots, where the critical temperature dependence on different coupling 
strengths can be easily seen. For the special case when both baths have the same temper¬ 
ature, we show that the critical temperature, above which the system becomes separable, 
satisfies /3cW ~ 2 (l -|- 4cT/3a;^) ^ for weak oscillator-bath coupling, cj being the oscillator 
natural frequency. It is consistent with the general expectation that l3cU} ~ 0{1) in the van¬ 
ishing inter-oscillator coupling a limit. In the opposite limit, when the oscillator-bath cou¬ 
pling is strong, correction terms with bath cutoff frequency dependence will show up. This 
is a noteworthy point in a lesser-explored regime, namely, one needs to be mindful of the 
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choice of the environment cut-off frequency in the treatment of strong system-environment 
coupling. 

A cautionary remark is in place here about entanglement measures used for quantum 
systems at finite temperature: Although the Peres-Horodecki-Simon (PHS) criterion is 
totally valid to identify the existence of entanglement in a quantum system, it does not serve 
as a quantifiable measure. We find from explicit calculations that at finite temperature 
it does not necessarily vary monotonically with the parameters in our system, namely, 
the temperature or coupling constants. One should exercise caution in using the PHS 
criterion for a physical understanding of thermal entanglement. In contrast negativity is 
a valid measure to quantify the dependence of quantum entanglement on these physical 
parameters. 

1.4 Differences from the common bath case 

To highlight the qualitative features in the behavior of the separability criterion it is useful 
to contrast the private bath case (Case Cl) studied here and the shared bath case (Case 
B) studied in Paper I. A more detailed description can be found in the last section: 

1. The initial Gaussian conditions will be irrelevant in the private bath case, but remain 
significant in the shared bath case, so the state of entanglement is sensitive to the 
initial conditions. 

2. At late times the entanglement measure for the private bath case is time-independent, 
but it continues to oscillate in time. 

3. The inter-oscillator coupling (a > 0) plays a more important role in the private bath 
case than in the shared bath case. 

4. In the private bath case, entanglement is easier to survive for stronger inter-oscillator 
coupling and weaker oscillator-bath coupling, but in the shared bath case, both factors 
seem to be overshadowed by the intrinsic quantum dynamics of the system which 
depends on the initial conditions of the oscillators. 

This paper is organized as follows: In Sec. 2, we briefly discuss the dynamics of the 
reduced system in the NESS configuration, and introduce the separability/entanglement 
criterion. In particular we pay attention to the covariance matrix, which constitutes the 
building blocks of the separability criterion. In Sec. 3 and 4, we highlight the calculations 
of the covariance matrix elements at high, zero and low temperature cases. We further 
examine the temperature dependence of the covariance matrix elements and the validity 
of the relevant approximations in Section 5. In Sec. 6 we investigate the separability 
criterion at different temperature regimes in detail and point out its non-monotonic be¬ 
havior. Because of this we adopt instead negativity as a valid measure of entanglement for 
quantitative analysis of quantum systems at finite temperature. We derive some relations 
between the critical temperature and various coupling constants. In Sec. 7, we then offer 
a more intuitive viewpoint to understand how all sorts of interactions can affect entangle¬ 
ment between oscillators. Finally we summarize our results and compare them with the 
case of the shared bath in Sec. 8. 
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2 The Model and the Covariance Matrix of the System 
2.1 The Model 


Consider two coupled harmonic oscillators of equal mass m and (bare) natural frequency 
Ub coupled to each other with strength a, each of which interacting with its own thermal 
bath with coupling constant e. We refer to the two oscillators together as the system, and 
the two baths together as the environment. This setup is a prototype used often for the 
investigation of nonequilibrium steady state (NESS), the existence of which is shown in a 
recent paper [1] (see also the references therein). In the Langevin equation approach the 
two oscillators’ amplitude Xi,X 2 satisfy the following equations of motion: 


Xi + 27X1 + w^Xi + 0 -X 2 
X2 + 27 x 2 +uj‘^X2 + crxi 




( 2 . 1 ) 

( 2 . 2 ) 


where 7 is the damping constant related to e by 7 = e^/(87rm), and uj is the renormalized 
frequency (wherein the correction from the interaction with the environment has been 
considered before), and ^ 1 , ^2 are the stochastic forces acting on Oscillators 1, 2 (© 1 , 2 ) 
respectively. Note they are not specified by hand but determined self-consistently. An 
overdot denotes taking the time derivative of a variable. The initial state of the oscillator 
is described by a Gaussian wavepacket and both oscillators are prepared in the same initial 
configuration. The two private baths (i?i, 2 ) are modeled by massless scalar fields at different 
temperatures /3“^. 

In the matrix notation, these two Langevin equations are condensed into one, namely. 


where 


X + 27X + n2 -X 


1 

m 




X = 




The solutions to this equation are given by, 


X{t) = Di(t) • x(0) + D 2 (t) • x(0) + — / ds D 2 (t - s) • ^(s). 

m Jo 


(2.3) 


(2.4) 


(2.5) 


where x(0), x(0) represent the initial configuration of the oscillators. The fundamental so¬ 
lution matrices Di, D 2 are a special set of homogeneous solutions to the Langevin equation 
(2.3), 


Di(0) = l, Di(0) = 0, D2(0) = 0, D2(0) = 1. (2.6) 


In particular, the Fourier transformation of 


(2.7) 
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is equal to 6{t) D 2 (t), that is, 


6{t) D2(t) 


dn 

-^21+02 _i2Ki ■ 


( 2 . 8 ) 


The function 0{t) is the unit-step function. Unless mentioned otherwise, we will not 
distinguish 9{t)Y)2{t) from D 2 (r) for all practical purposes, and denote (—— 
i2 acI) ^ by D 2 (k). 

The force term ^i{t) is a stochastic c-number with the statistical properties 


(^W) = o, 


=e^ Gnit-t') = e^ 


0 Gjf(t-t')J’ 


(2.9) 


where G^^{t — t') is the Hadamard function of the bath scalar field, associated with the ith 
oscillator [1]. This stochastic force in essence represents the quantum fluctuations of the 
private bath at a finite temperature. 


2.2 Entanglement Measures 

For continuous-variable systems, the entanglement measure based on the density matrix is 
not conveniently calculable because the density matrix in this case is infinite-dimensional. 
However, it has been shown [18] that in the case of continuous Gaussian variables, the 
Peres-Horodecki separability criterion [16, 17] can be reformulated in terms of the covari¬ 
ance matrix of the bipartite system, 

C+ = det A det B - Tr{ A • J • C • J • B • J • • J} + (det C + 

— - (det A -|-det B) > 0 , (2.10) 


with 


J = 




Here the matrices A, B, C are the block matrices in the covariance matrix V, 


V = 



( 2 . 11 ) 


while the covariance matrix V itself is defined by the canonical variables of the two sub¬ 
systems 


V = - Tr 
2 


p{R,R^}J =-({R,R'^}), 


( 2 . 12 ) 


where p is the density matrix of the state we are interested in. We have assumed (R) = 
0. The column matrix R takes the form R^ = (xi,pi,X 2 )P 2 ); and pi is the canonical 
momentum conjugate to Xi associated with the subsystem i. The angular brackets denote 
taking the quantum expectation value. In our case, once we have the covariance matrix for 
the coupled harmonic oscillators in the NESS configuration, we may construct C+ according 
to (2.10). A negative value of C+ thus implies the existence of quantum entanglement. 
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Although (2.10) constitutes only the second moments of the canonical variables, it 
offers a complete description of the Gaussian system since for a Gaussian system, all higher 
moments can be expressed in terms of the second moments. Oftentimes it is instructive 
to write the Peres-Horodecki separability criterion in terms of the symplectic eigenvalues 
of the partially transposed covariance matrix. Let r]^ stand for the symplectic eigenvalues 
of the partial transposition of V. Without loss of generality we assume rjy is greater 
than 7y<. In fact they can be found by solving the eigenvalues of the matrix ifl ■ with 
fi = ©LiJ- The resulting eigenvalues will appear in pairs by the form ±r/>, ±^<, so the 
symplectic eigenvalues of are given by the absolute value of the eigenvalues of z O • 
When we write into the Williamson’s form, the separability criterion + zfi/2 > 0 
becomes 


/ v> 

0 

0 

0 \ 


( ° 

+1 

0 

0 \ 

0 

r}> 

0 

0 

i 

+ 2 

-1 

0 

0 

0 

0 

0 

V< 

0 

0 

0 

0 

+1 

Vo 

0 

0 

V<) 


VO 

0 

-1 

0 / 


0, ^ (2.13) 


When ry< < 1/2, entanglement occurs. Notice that r/> is assumed to be larger than r/<, so 
r 7 > is always greater than 1/2. We observe that although a violation of the Peres-Horodecki- 
Simon separability criterion signals the existence of entanglement, it is not a good measure 
for a quantitative description of entanglement, in that the criterion includes a unwelcome 
factor (?]> — 1/2), which does not affect the identification of the existence of entanglement, 
it messes up the correct evaluation of entanglement. This can be understood if we examine 
the behavior of the symplectic eigenvalues about r^< ~ 1/2. For definiteness, we assume 
that the symplectic eigenvalues are similar monotonic functions of the parameters of the 
entangled system. We can easily see that if t/> changes too fast in the vicinity of r^< = 1/2, 
the product (r/^ — 1/4) (r/^ — 1/4) will not be monotonic there. 

As is perhaps better known, a simple calculable measure of entanglement which also 
provides quantifiable accuracy is negativity [20], denoted by M or its logarithm (strictly 
speaking logarithmic negativity is not merely the logarithm of negativity, although it is 
related to) [21], the logarithmic negativity E/,/. For the Gaussian states under study they 
can be respectively defined by 

AA(p) = max{0, , Ex{p) = max{0,-In 2r]<} , (2.14) 

zil< 

in terms of the symplectic eigenvalue rj^ of the partially transposed covariance matrix. 
When < 1/2, the Gaussian state p is entangled and both measures take nonzero values 
between O"*" to -|-oo. In addition, the logarithmic negativity has a convenient feature of 
being additive. 

Gomparing the negativity (2.14) with the Peres-Horodecki-Simon criterion (2.13), we 
observe that they are all based on the smaller symplectic eigenvalue ?/< of the partially 
transposed covariance matrix V^*, so they will give the same prediction on the existence 
of entanglement. However, the separability criterion carries an additional undesired factor 
(r/> — 1/2), which may inadvertently scale (r/< — 1/2). Thus the separability criterion is 
not suitable for quantifying entanglement. 
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Finally, we remark on a subtlety of the entanglement measure. It has been pointed 
out [22-24] that different measures may give different ordering of density operators with 
respect to the amount of entanglement. To be more specific, given two density matrix pi 
and p 2 , we can have Ei{pi) < Ei{p 2 ) for one entanglement measure, while E 2 {pi) > E 2 {p 2 ) 
for another. In particular, negativity and Gaussian entanglement of formation, the latter 
forming an upper bound to the true entanglement of formation, have been found to be 
inequivalent for asymmetric Gaussian states [24]. For symmetric states, the predictions 
from both measures coincide. 

The next few sections will be dedicated to the calculation of elements of the covariance 
matrix. 

2.3 Elements of the Covariance Matrix 

We use (2.5) to find the elements of the covariance matrix V. Assume that the initial state 
of each oscillator is depicted by a Gaussian wave packet of the same shape, at rest initially 

at the bottom of the harmonic potential associated with each oscillator, such that 

{x^m = (MO)) = 0, {{MO),MO)}) = 0, (2.15) 

({w(0),Xi(0)}) = (x?(0)) 6ij , {{Pii0),Pj{0)}) = {PiiO)) Sij , (2.16) 

with Pi = mxi- Thus these two oscillators are initially in a separable state. From the 
solutions (2.5) one can identify the role of the interaction, either between the oscillators or 
between the oscillator and its private bath, in creating or sustaining the quantum entan¬ 
glement in the system. 

To calculate the elements of the covariance matrix V one can show, for example, that 
!({»(*), xj(t)}) = D;‘(t)Df (*) {xlm + Td?(«)d?(() {plm 

+ ^ [ dsds'Bf{t-s)Bi\t-/)G’lf{s-s'). (2.17) 

Jq 

When the dynamics of the system evolves into relaxation as f —?• oo, the first two terms 
on the righthand side will be exponentially small if the coupling constant between the 
oscillator and the bath is not vanishing. Thus at late time ({xj(t), Xj(^)})/2 simplifies to 

lim -{{xiit),Xj{t)}) = W 2 / (s') Gh{s - s') 

t^oo z m j—QQ 

2 foo _ 

= («) cjf(«), (2.18) 

J —00 

where we have used the fact that D 2 (r) = 0 if r < 0. Since the Fourier transform of D 2 (s) 
is defined by 

^ + -ilKl ’ 

we use the property D 2 (—k) = D 2 (k) to arrive at (2.18). 


(2.19) 



At this point, let US look at a more specific example: the element Vii(t) = Xi(^)})/2 

(Xi(t))- At late time it takes on the value Vn, 


„2 poo j 

Vn= limyn(t) = ^ / 

t^oo J_^ 2 tt 


|D^i(«)pG}i(K) + |D^2(^)pGi(K) 


( 2 . 20 ) 


and 



{k? - 

[(k^ — oj\Y + 47^k^] \{k^ — ’ 


[(k^ — + 472K2j _ ^2 ^2 _|_ 4^2^2j ’ 


( 2 . 21 ) 

( 2 . 22 ) 


with 7 = e^/(87rm). The frequencies ± u are the oscillating frequencies of the 

normal modes, which can be constructed from the superpositions of (2.1) and (2.2). The 
Fourier transformation of the Hadamard function takes the form 



K 

dvr 



K K 

- -1-e 

dvr 27 r 

1 


-/3fcK 


27r/3fc ’ 


/IfcK > 1 , 
l3kK <C 1. 


(2.23) 


The term K/dvr represents the vacuum zero-point contribution. The off-diagonal terms of 
Gh are zero because both private baths are not correlated. 

From the late-time value Vn of the amplitude uncertainty of Oi, we observe the fol¬ 
lowing distinct features: (1) it approaches a constant independent of time, (2) its integral 
expression (2.20) takes a form similar to the Landauer formula, where |D 2 ^(k)P plays a 
role of the transmission coefficient, and (3) it depends on both thermal baths even though 
Oi does not have a direct contact with 32- The last property would not be unexpected 
because the coupling between the oscillators will bring in correlations between Oi and B 2 , 
and vice versa, between O 2 and Bi. In fact, these features hold quite generally for the all 
elements of the covariance matrix in the current NESS configuration. 

The definition of the covariance matrix V and the expressions for its elements, and 
their corresponding values at late times are derived in Appendix A-D. 

In the next sections we will explicitly evaluate the elements of the covariance matrix 
for three situations: (1) high temperature limit, (2) zero temperature case and (3) low 
temperature regime. 


3 The Covariance Matrix at High Temperatures 


We consider the high temperature limit jSoj <C 1 of the elements of the covariance matrix. 
In this limit, the Hadamard function of the bath (2.23) is approximately given by 


G“(k) c. 


I 

2TTf3k ’ 


(3.1) 


We see that its vacuum contribution is relatively negligible, and can be neglected for most 
cases. However, extra discretion is advised for the evaluation of the momentum uncertainty 
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where the vacuum contribution of the bath can be significant when the coupling between 
the oscillator and the bath is sufficiently strong. Thus the result can depend on the cutoff 
scale of the environment field) (see, e.g., [12]). 

Here we merely highlight the calculation for the element Vu at late time. To obtain 
the high temperature limit of Vii(oo), that is, Vn, essentially we evaluate the following 
two integrals 

vr / 47^ \ 

47 4^272 -g o- 2 y ’ 

(3.2) 

TTCT^ /I 1 \ 

4^27 ^ 1^4 _ ^2 4^2^2 _|_ ^2 J • 

(3.3) 


/ OO 

dn 

-OO 


(^2 — (^2^2 _|_ 4^2^2 


[(^2 — a;^)2 + 472K2j [(^2 — (^2 ^2 _j_ 4-y2^2j 


h = i dn 


a 


[(^2 — a;^)2 + 472 ^ 2 ] [(^2 — a;2 )2 + 4 ^' 


2^2 


In terms of I\ and /2, we see from ( 2 . 20 ) that the high temperature limit of the element 
Vn at late time is given by 


27 

7rm [/3i P2_ 

1 f 8 w"^7^ + - 47 ^ 0 -^ 1 cr^(w^+47^) 1 1 

2 m \ (w^ — + a'^) fdi (w^ — cj2)(4a;2-)/2 _|_ (j2j ^2}' 

Here we would like to point out that when the mutual interaction a is large, in particular 
when a ^ uj, the fluctuations of the oscillator grow significantly. This will be traced back 
to the small values of a;_. We will come back to this feature in due course. 

Derivations of the high temperature forms of V13, Vu, V22, V24 are given in Appendix 
B. Nonetheless for the following discussions we will bring forward the results for V22 and 
V13 here. When both private baths have the same temperature /d~^, we have from ( 5 . 1 ), 
( 5 . 3 ) and ( 5 . 5 ) 


Vn = (x?(oo)) = 


Vn 

V 22 

Vl3 


(x?(oo)) 

(Vi(oo)) 


1 

m/3 — 0-2 ’ 

m 


(xi(oo)x2(oo)) = 


1 a 

m/3 — cj2 ’ 


(3.5) 

(3.6) 

(3.7) 


in the weak oscillator-bath coupling limit. This implies that the average harmonic potential 
energy of Oscillator 1 (Oi) will be 


Esi 


m 

~2 


‘^^(x?(oo)) 


1 

^ a;4 - 0-2 ■ 


(3.8) 


It is a bit off from the value 1 / 2/3 one would expect from the equipartition theorem for a 
free harmonic oscillator in the high temperature limit. This difference will disappear when 
the mutual coupling cr between the two oscillators are turned off. 
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Eq. (3.6) on the other hand tells us the corresponding average value of the kinetic 
energy in the high temperature limit, 

_ V22 _ 1 
2m 2/3 

is the same as the value obtained from the classical equipartition theorem. We observe that 
in the high temperature limit the mean kinetic energy is not equal to the mean harmonic 
potential energy in general, and the sum of the kinetic energy and the harmonic potential 
energy is not equal to kT: 

Efc,+E,, = —+ —(3.10) 

Let us compare this with the average total energy of a free harmonic oscillator in a closed 
system, 


{H) =Ek + Es 


00 


Y,Ene-^^' 

n=0 


n=0 


dlnZ 1 


00 


Z = Y^ e-^^' 

n=0 


(3.11) 


in the high temperature limit and = (n + ^)uj- 

The deviation can be accounted for by the fact that some portion of the total energy 
of both oscillators is stored in the mutual interaction between Oi and O 2 . Accordingly the 
missing piece should come from the expectation value of ma xiX 2 - Its contribution to the 
mechanical energy is 

E^ = lim ma {xi{t)x 2 {t)) = maVis = —2 ’ 

t^oo p (jj^ — a^ 


when = P = p 2 - Including this contribution we see the total energy for the two-oscillator 
system in the high-temperature limit becomes 


E — Ek^ + Esi + Ek^ + Es2 ~ ~p ~p ^4 _ fj2 


1 


a 


Pio^-a^ p 


(3.13) 


which is that obtained by the classical equipartition theorem for two coupled linear oscil¬ 
lators. This also serves as a consistency check of our calculation. 

Finally we comment on two issues. First, weak oscillator-bath coupling enables us 
to ignore the cutoff-dependent effect from the bath. This may not be true in the strong 
coupling case. Second, despite the resemblance of (3.13) with (3.11), they are quite different 
in the physical context. The former is considered in the context of open systems while the 
latter under the assumption of a closed system. It has been shown [13] that both results 
can be equivalent only in the limit of vanishingly weak oscillator-bath coupling. 
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4 The Covariance Matrix at Zero and Low Temperatures 

Here we evaluate the vacuum contributions and the low temperature correction of the 
covariance matrix elements. Due to the zero-point fluctuations of all bath modes, the 
vacuum contributions of some covariance matrix elements can be divergent. Suitable cutoffs 
need be introduced to regularize them, with due consideration of the particulars of the bath 
the system interacts with. 

Let us examine, for example, Vn = Vii(oo) and work out its zero and low temperature 
expressions. 


4.1 Vii at zero temperature 

The vacuum contribution of G^{k) is 

G“(k) = sgn(K) ^ , 

so we need the following two integrals to evaluate the vacuum contribution of Vn, 

-I- 47^k^] 


Ji = 
J2 = 


L 


dK 

dn 


lo [(^2 - )2 + 472^2] [(^2 - w2 )2 + 472^2] ' 

The sum of Ji and J 2 can be expressed as 


TT 




/(n+) /(!i_) 


n+ 

where the dimensionless function f{z) is defined by 

27Z 


r >2 2 2 

S — 7 


f{z) = 1 - 1 — cot 


TT z"^ — 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 


It is clear that rZ± are the resonance frequencies of the two normal modes. Therefore from 
(2.20), the zero-temperature (vacuum) contribution of Vn is given by 


, ,(0) _ ^ 

Kii — 

Trm 


Jl + J 2 

1 


8m 


/(H+) ^ fin.) 


Hi 


n. 


(4.6) 


We observe that the vacuum contribution can be clearly separated into decoupled com¬ 
ponents of two normal modes, with oscillating frequency oj± respectively. This is another 
general feature of this system. 

The zero-temperature expressions for V 13 , V 14 , V 22 , V 24 are given in Appendix C. 


4.2 Vn at low temperature /3a; /§> 1 

The low temperature corrections to the covariance matrix elements basically result from 
the corresponding correction in the Hadamard function, 

00 

G“(K)c^vac. + |^ (4.7) 

n=l 
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because the fundamental solution matrix Di ^2 does not depend on temperature. This is 
a consequence of the fact that the retarded Green’s function of the scalar field, which 
accounts for dissipation in the Langevin equation, is state-independent. 

As is seen from (2.20), we need the following two integrals to evaluate the low temper¬ 
ature correction of Vn, 


Ki=2 
K2 = 2 


f 

f 


dn 


dn 


-|-47^k^] e 20 ;“^ 1 

[(k^ — a;^)2 -|- w+wi 

2 (t 2 1 


-|- 472^2] [(k^ — uj'^Y + 47 ^k^] w+wi fJ'^ 
We then have the low temperature correction to Vn given by 

1 


V 


(/3) 

11 


27 


47 

'iO^ 


rvj - 

K 1 + K 2 

= - 

nO 

+ 

Trm 


nm 




A, A 


+ 0(r'), (4.8) 
+ 0(^-3). (4.9) 


(4.10) 


However this is merely the contribution from the first term in the summation of all finite 
temperature corrections in (4.7). Since the remaining terms (n > 1) will have algebraically 
comparable contributions, we have to take them into consideration. We note that the 
leading term in has a temperature dependence /3~^ in the low temperature limit. 
Thus we expect that the leading contribution for the remainder of the finite temperature 
corrections in (4.7) should be proportional to Their overall contributions will 

introduce an addition factor 


00 ^ 

n^/3^ 

n=l 


TT^ 1 
~6^’ 


(4.11) 


to (4.10). Therefore after taking this into account, we obtain the low temperature correction 
to Vii as follow: 


V 


(4) 

11 


2 Try 

3m 


7 ^ 


4 4 

cjIcjI 


(4.12) 


The low temperature expression of Vn is then given by the sum of (4.6) and (4.12), 


vu = vi;>+vif = i 


7(^+) , /(^-) 




3m [jJl ^ Pi 


2tt^ 


4 4 
wlcu 




(4.13) 


We leave the derivations of the zero and the low temperature expressions for V 13 , V 14 , V 22 , 
V 24 in Appendix D. 


5 Temperature Dependence of the Covariance Matrix 

Because elements of the covariance matrix at finite temperature may prove useful for more 
general purposes, we collect their expressions for both high and low temperatures for the 
system at late times when it reaches a NESS, the existence of which for this setup is shown 
in our earlier paper [1]. 

Here we summarize the temperature dependence of the elements of the covariance 
matrix. 
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1 - 1^11 = 


Vii = 


f{n+) ^ f{Q.) 




0 _ 


+ 


27r7 io'^ 

3m [pf Pi 


(u'^ — 0 - 2)2 


1 

8m 
1 

2m [(o;^ — o- 2 )( 4 u ;272 + 0 - 2 )/?! ' — o- 2 )( 4 a ;272 + cr2)/32 


8w'^ 7^ + w^o-^ - 47^0-^ 1 o-^(w^ + 47^) 


/ 372 W > 1 , 

, /3i,2W <C 1 . 

(5.1) 


2 . V12 = -({7i(oo),pi(cx))}): 


V12 = 0 . 


(5.2) 


3 - Vi3 = -({71(00),72(00)}): 


1 



2777 Ol^fJ 

r 1 11 

I 

8m 


0 - 

3m ( 07 ^^ — 0-2)2 



Vl 3 = 


1 a 

2m — 0-2 


1 1 

A 


/3i,2W >• 1 


P 12 OJ <C 1 . 


(5.3) 


4 . Vi4 = -({71(00),p2(oo)}): 


Vl 4 = 0 - < 


Svr^ 7 ^ 0 - 1 1 

1^(u;4-ct2)2 

70- 


4^272 + 0-2 


1 1 

/3i /32 


/3i,2<^ ^ 1 , 
Pi,2^ <C 1 . 


(5.4) 


5. V 22 = -({pi(oo),pi(oo)}): 


m 

271 o;^ — 0-2 8 


rn'y 

V 22 = — In ^ 


— 7 ^ — 7 ^ 


+ 


471^ 




m7 




O" 

TTT A ^TT 


15 (u;4_^2)2 L^4 ^4J ’ 


mj 

TT 


^0(/?,A-l) ln(/3,A) 


m 


8a;^7^ + 0-2 1 


+ 


/3i,2W>1, (5.5) 

2 1 - 


a 


407272 + (j2 407272 + (j2 ^2 

/3i,2W <C 1. (5.6) 


Mathematically speaking, the inclusion of the unit-step function 0(/3A — 1) is to ac¬ 
count for the vacuum contribution of the bath modes in the case /3A > 1, because 
■when Pk > 1, the Hadamard function G^{n) takes the low-temperature form as 
shown in (2.23). On the other hand when pA < 1, the high-temperature approxima¬ 
tion of G^{k) is entirely valid up to the cutoff scale, the cutoff-dependent term being 
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subdominant. However on physical grounds, since the cutoff scale by construction is 
the highest energy scale compatible with the model, the thermal excitation energy 
thus must be smaller than the cutoff scale. It then implies that even in the high 
temperature limit, we still have oj <C j3~^ < A. 




Here we show the high/low temperature approximations of V 22 with a numerical 
calculation. In particular we explicitly highlight the role of the vacuum contribution, 
that is, the cutoff dependent terms, even in the high temperature approximation for 
strong oscillator-bath interaction. The pink curve in the plot on the right shows 
that if the vacuum contribution of the bath is not taken into account, the analytical 
high-temperature approximation will be way off from the numerical result (the purple 
curve) in the region /3a; ~ 0{1). On the other hand, the red curve, which includes 
the vaccuum contribution, fits nicely with the numerical result. The parameters are 
chosen to be 7 = 0 . 2 , a = 18, a; = 5, and A = 1000. The plot on the left is drawn 
for weak oscillator-bath interaction 7 = 0.2, i.e. 7 / 0 ; <C 1. The cutoff-dependence is 
seen as dispensable. 

6 . V24 = ^({pi(oo),p2(oo)}): 


,, m'y , 0 ;+ m 

V24 =-^ In — 


TT 


U!- 


O?. — 7 ^ 0^ — 7^ 


Q: 


47r^ 


muj ja 


15 (a;4-cT2)2 [pf ’ 


1 1 

TTT T WT 


ma 


/3i + /32 


/3i,2a; > 1, (5.7) 

/3i,2a;<Cl. (5.8) 


In this case, since the leading contribution of the high-temperature approximation 
vanishes, we have to include the next-order term. 

1 


7 . V23 = -({x2(oo),pi(oo)}): it is equal to -V14, so 


V23 — — Vl4 — 0 < 




7 ^cr 


1 


1 


15 (a;4-cT2)2 Ipf 

7(T 


4a;272 -|- (j2 


1 1 

/3l f^2 


/3i,2w 1 , 

'C 1. 


(5.9) 
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8- '1^34 = -({X2(oo),P2(oo)}): 


V 34 = 0. 


(5.10) 


9. V 33 = -({x 2 (oo),X 2 (oo)}): it is similar to Vn except that we swap /3i and (32, 


V33= < 


1 

8 m 

1 


f{n+) ^ f{Q.) 






+ 


27r7 

3m 


ui^ a 




— 47^(7^ 1 


2^2 


+ 47 ^) 


2 m [(o;^ — + cj^) /S 2 — (T 2 )( 4 ti ;272 + u^) /3i J 


/3i,2W » 1 , 

/3i,2W <C 1 . 
(5.11) 


10 . V44 = 2(i^’2(oo),p2(oo)}): it is similar to V22, 


rri'y 

V44 = w— In ■ 


27r uj^ — a‘ 


m 

+ TT 


111 - 7 ^ , 111-72 


Hi 


/( 11 +) + 


H_ 


m-] 


+ 


47r^ 


mj 


15 (a;4-a2)2 [/34 ’ 


1 ! 1 ! 

4~ 


/?i, 2 a;>l, (5.12) 


m7 

TT 


j;0(/3,A-l) ln(/3,A) 
i=i 


m 

+ — 


8lo^Y + (^ 1 


+ 


a 


2 [ 4 a;V + f ^^/52 4^272 + / 3 iJ ’ 

/3i,2W <1. (5.13) 


Some comments are in place here: Both oscillator are initially prepared in a state of non¬ 
overlapping Gaussian wavepackets with the same width <;■. As they come into interaction 
with their own private baths, the evolution of each individual oscillator is then driven by its 
bath and the other oscillator it is directly coupled with. Due to the dissipative self-force on 
the oscillator arising from its interaction with its own bath, the intrinsic information of the 
initial state is dispersed away exponentially fast as the system evolves in time. In the end 
when t —)• 00, the values of the dynamical variables of the oscillator are determined by its 
private bath and by the other oscillator. We want to bring up this point because even when 
the oscillator-bath coupling constant 7 approaches zero, not all of the asymptotic values 
of the covariant matrix elements are zero. In this limit their values are independent of the 
parameter characterizing the initial state, so they are not related to the intrinsic evolution 
that begins with the initial configuration. Instead they are the induced components as a 
consequence of the interaction between the oscillator and the bath. In other words, the 
results of the covariance matrix in the vanishing 7 limit should be understood by the 
limiting procedures of taking t —)■ 00 first and then taking 7 — 7 - 0 . 

This is a good point to discuss in more details in what is meant by the high/low 
temperature approximations. We only cover the generic situation and discard some extreme 
cases such as uj-, H_, cr —)• 0 , so we assume that era, and Hj- are about the same order 
of magnitude as the oscillating frequency uj. The cutoff frequency is assumed to be much 
larger than uo, i.e., uj. The magnitude of the parameters 7 and are not restricted 
as long as H_ remains well-defined. We use Vn and V 22 as illustrating examples. 
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1 . Vii: as far as the order of magnitude is concerned, we see 


vacumm: 

low temp: 
high temp: 


(0) ^ ^ r/^ /(^-) i ^ j_ 

8m mo; ’ 

(^) ^ 2^\u^ ^1 1 _ 1 7 1 

3m /3f /3| (o;^ — fx^)^ moj oj {l3u:Y ’ 

_ 1 8a;^7^ + — 47 ^ cr^ 1 (T^(w^+47^) 1 

2m (o;^ — + cj^) /3i (o;^ — (T 2 )( 4 ti ;272 + (52 

1 1 

rs^- 

muo l3oj 


Roughly speaking, the high temperature limit refers to the case /Jcu <C 1; on the other 
hand a consistent low temperature correction requires 


7 1 

W (/?W)2 


< 1 , 


which can be weaker than the naive low temperature limit /3a; ^ 1, especially in the 
weak coupling limit 7 / 0 ; <C 1. It implies that in the weak oscillator-bath coupling 
limit, the low temperature correction has a much wider range of validity. In the 
strong coupling regime 7 ^ 0 ;, the low temperature correction is remains fully valid 
for the /3a; » 1 . 


2. V22: 


vacuum: 




Il 2 _^2 


7 A 
muj — in¬ 


ert) 47r^ 

low temp: V22 = 


Itt'^ m7 lo;^ cj^ ] 7 1 

(a;4 - n 2)2 ^ ^ ’ 


1 . , ,,(/ 3 ) fn 80^272-g 0-2 1 0-2 1 1 

high temp: Voo = — . o o - 9 TT + . 00 -oTT ^ muj — . 

22 2 [ 4 a ;272 + 0-2/ 3 i 40^272 + 0-2/32J I 3 u } 

Here, additional subtlety arises due to the presence of the cutoff frequency A. The 
importance of the cutoff-dependent term relies on how the factor 


is compared with unity. In the weak coupling limit, the cutoff dependent term is 
negligible, so we can safely ignore it unless the cutoff frequency is extremely high, 
such as 

A ~ O{coe^). 


In the strong coupling regime 7 < a;, we see that the cutoff-dependent term still 
has a comparable magnitude relative to the high temperature approximation in the 
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interval of the high-to-low temperature transition j3ui ~ 0{1). This interval has a 
special significance because, as we will see later, this is the region where thermal 
entanglement may disappear in the nonequilibrium steady state configuration. 

Thus at this point, generically speaking, the high-temperature limit refers to j3uj <C 1 
while the low-temperature limit refers to /3a; ^ 1. For weak oscillator-bath coupling, the 
low temperature correction has a wider range of validity than is implied by /3a; S> 1 due 
to the additional factor 7 / 0 ; in the corresponding expression. In addition, the cutoff is 
completely negligible in normal circumstances. By contrast, in the strong coupling regime, 
the cutoff-dependent contribution enters in determining the critical temperature of thermal 
entanglement. 


6 Entanglement of System in Noneqnilibrium Steady State 
6.1 Late Time Behavior of the Covariance Matrix 


At late time when the system reaches the steady state, the covariance matrix takes the 
form 

/Vii 0 Vi 3 Vi 4 \ 

0 V22 ~Vl4 V24 

Vl 3 —Vi 4 V33 0 

\Vl 4 V24 0 V 44 y 


with 


V = 


/Vii V12 Vi 3 Vi 4 \ 

V21 V22 V23 V24 
V31 V32 V33 V34 
\^V 41 V42 V43 V 44 y 


A C 
B 


A = 


C = 



Vl 3 Vl 4 
—Vl 4 V24, 


B = 


J = 



( 6 . 1 ) 


( 6 . 2 ) 

(6.3) 


The determinants of the matrices A, B are related to the generalized uncertainty relation 
for each single subsystem, which also takes into account the correlation between canonical 
variables. The matrix C contains the cross-correlation among canonical variables between 
two subsystems. 

As is briefly discussed in Sec. 2 . 2 , the knowledge of the covariance matrix enables us to 
use the Peres-Horodecki-Simon separability criterion to determine the quantum entangle¬ 
ment. In fact the separability criterion can be combined with the generalized uncertainty 
relation to form an unihed statement 

C± = det A det B - TV{ A • J • C • J • B • J • C'^ • J} + (det C ± ^ (det A + det B) > 0, 

(6.4) 


The expression containing the — sign represents the uncertainty relation while that with 
the -|- sign represents the separability criterion. We immediately see that in the current case 
det A and det B are always positive dehnite by construction. In addition, the expression 
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Figure 1. The separability criterion C_|_ is plotted against the oscillator-bath coupling 7 and the 
mutual coupling strength a between the oscillators. The black curve demarcates the separate 
(C+ > 0) and the entangled (C+ < 0) regions. The choices for the parameters are w = 5 and 
A = 10000. 


Tr{A J C J B J J}, once written explicitly in term of the covariance matrix 

elements, 

Tr{ A • J • C • J • B • J • • J} = V22V44Vf3 + V11V22VI4 + (V11V44 + V22V33) > 0 , ( 6 . 5 ) 


is found to be always positive. Thus when we rewrite ( 6 . 4 ) as, 


c± = 


1 


det A —- det B —7 + det C ± -- 


)+(deter 1)' - Tr{ 


TV<^ A J C J B J 


T + ie 

(6.6) 


we immediately recognize that actually contains two positive but competing compo¬ 
nents. This makes it difficult to determine the sign of Ci- However, we can use the 
following argument: Suppose that the uncertainty relation C- ^ 0 always holds. Since 


C+ = C- + det C , ( 6 . 7 ) 

the condition ^_ < 0 implies that det C must be negative. Therefore the appearance of 
negative values of det C may help to signify the existence of entanglement. The sign of 
det C is less clear, 

det C = V13V24 + Vu , ( 6 . 8 ) 

depending on how negative V13V24 can be allowed. Although this is not a sufficient condi¬ 
tion, it highlights the role of cross-correlations in entanglement. 

Before we proceed to evaluate C+, we observe that among the elements of the covariance 
matrix, two of them, V22 and V44, have dependence on the cutoff frequency A, which is 
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the highest energy scale that is consistent with the theory. Thus we expect that C±) and 
in particular, the separability criterion, will depend on the cutoff scale. Since the cutoff- 
dependent term always has the form 7 In A, where 7 is the system-environment coupling 
constant, it implies that this cutoff dependence will be suppressed in the weak coupling 
limit. However, when the system interacts strongly with the environment in the sense that 
7 / 0 ; is close to unity, the contribution from the factor 7 In A can be significant, and can 
make the separability criterion ambiguous. 

Likewise in terms of the matrices A, B, and C we can construct the symplectic eigen¬ 
values rj^ of the partial transpose of the covariance matrix V [25], 

1 
2 

, ( 6 . 


= 


det A -|- det B 


- - det C ± 


’ det A -|- det B 


-detC -detv 


where alternatively det V can be written as det A det B -|- (det C)^ — Tr{ A • J • C • J • 
B J J}. This enables us to calculate the quantitative entanglement measures like 

negativity or logarithmic negativity for the Gaussian state. 

In the sections that follow, we will refer to the special case when both thermal reser¬ 
voirs have the same temperature. In this case, the Gaussian state becomes symmetric, so 
(logarithmic) negativity will give an unambiguous ordering of density matrices, in com¬ 
parison with other quantitative entanglement measures. Since we have A = B, and the 
matrix C becomes diagonal, the symplectic eigenvalues 77 ^ takes a particularly neat form 

( 6 . 10 ) 


= 


(Vii T Via) (V22 ± V24) 


with 

Vii = (A?(i)), V22 = (p?(i)), Vi3 = ^ ({xi(t),X2(t)}), V 24 = ^ ({pi(t),P2(i)}) • 


We readily see that 

Vii =F Vi3 = ^ ({xi,Xi =F X 2 }), (6.11) 

V22 ± V24 = - ± P2}), ( 6 - 12 ) 

are associated with the dynamics of the normal modes of the joint system. This elicits a 
transparent connection between the entanglement behavior and the underlying dynamics. 


6.2 Entanglement Behavior 

As stated earlier the Peres-Simon-Horodecki criterion can be used to identity the existence 
of entanglement of the Gaussian states, but it may be inadequate to provide a quantitative 
description of entanglement, in particular, for a quantum system at finite temperature. 
We will show later that it does not vary monotonically with temperature and coupling 
constants. This discrepancy comes from an additional factor in the criterion. It has no 
effect on the identification of entanglement but it will give an unwarranted bias on the 
values, rendering it inappropriate for quantifying entanglement. While the separability 
criterion can be used for the system under study at zero temperature, we need a different 
measure to quantify finite-temperature entanglement, namely, negativity. 
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c+ 



Figure 2. The separability criterion C+ is plotted against the mutual coupling strength cr between 
the oscillators at zero temperature. Larger values of the damping constant 7 will move the curve 
upwards and make the entanglement between the two oscillators harder to sustain. The oscillating 
frequency w and the cutoff frequency A are chosen to be 5 and 10000, respectively 


6.2.1 zero temperature 

We first examine the separability criterion at zero temperature. The whole expression 
for can be exactly found but it is tremendously large. Here we present the leading terms 
in the weak oscillator-bath coupling limit, i.e., 7 < oj± is the smallest parameter at hand, 


{UJ+-U}-Y {U)+-U]-) 

IGcu+cj- 

,.,2 

+ 


2 2 ^ I 

w, — 0J_ + In 

0J_ 


7 


327r2a;^a;2 


ITT^ (a;+ — a;_) ^ + 4a;+w_ + ) — 32a;^a;^ 


-b IGcj+W- 


UJ- 




A A 


A 


A 


2 I 

uj_ In — 
OJ- 


+ 0(7lnA)3. 


(6.13) 


Note that it is not sufficient to expand to first order in 7 because they all depend on 
(a;+ — oj-). This factor will make the first-order expansion of ()+ vanish when cr —0 no 
matter what value 7 has. In fact ((7 has a finite value when 7 7^ 0 , so we have to include 
terms which are at least of second order in 7. 

In addition, as far as the leading contribution is concerned, we see that (a;+ — is 
always positive, so C+ is negative for all nonzero mutual coupling strength cr between the 
two oscillators. This implies that the oscillators will become entangled aympotically once 
they are coupled. On the other hand, when we consider contributions due to the finite 
value of the damping constant 7, we find that, in particular in the limit cr —?• 0, we have 
a;+ —)• uj- and 

72 / A \ 2 

lim C+ = ( In- 1 ) >0. (6-14) 
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Figure 3. The separability criterion C_|_ plotted against the cutoff scale A. The oscillating frequency 
ui and the inter-oscillator coupling a are chosen to be 5 and 21, respectively 


The separability criterion is positive for cr = 0 when 'j ^ 0. With increasing a, the 
value of C-i- gradually falls below zero at some critical value of CTc. Thus we see that the 
curve of the separability criterion will move upwards with larger values of the damping 
constant 7 , that is, with stronger interaction between the oscillator and its private bath. 
Furthermore, it also indicates that these oscillators are not always entangled, and they can 
be separable for some choices of 7 and a. For a specific value of 7 , the mutual coupling 
strength a must be greater than the critical value to render both oscillators entangled. In 
other words, the bonding between two oscillators has to be strong enough to overcome the 
incoherent disturbance from their respective baths in order to maintain their entanglement. 
The larger the values of a the more stable the mutual entanglement is. Therefore we see 
that the oscillator-bath interaction and the coupling between oscillators play competing 
roles in sustaining the entanglement. 

We now derive a relation between the critical values of different couplings. For the 
case of a small damping constant 7 , the critical value Uc can be obtained by solving (6.13) 
with = 0 , yielding 


TTfj 1 vrcr^ 1 

4wln^-l 4^3 (In A _ 1)2 ^ 


(6.15) 


Inverting it leads to 
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r A 1 
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(6.16) 


Taking Fig 2 as an example, setting a; = 6 , 7 = 0.135 and A = 10000 in (6.16) gives = 
5.868. It is nicely consistent with the intersection of the green curve with the horizontal 
axis in Fig 2 . 

Therefore in the weak coupling regime where ^ In ^ is small but not vanishing, we hnd 
that if the two oscillators are initially prepared in a disentangled state, they can become 
entangled for sufficiently strong direct mutual coupling between the oscillators. Otherwise, 
they remain asymptotically in a separable state when the mutual coupling is weak. Finally, 
we add some remarks on the cutoff dependence of the separability criterion. From (6.16), 
we see that the dependence on the cutoff always occurs as long as 7 7 ^ 0. This implies 
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Figure 4. The separate criterion C_|_ is plotted with respect to the temperatures of two private 
baths. The black curve C+ = 0 divides the separable state (C+ > 0, pink shade) from the entangled 
state (C+ < 0, green shade). Essentially the curve traces along the region fioj = 0(1). The 
oscillating frequency uj and the cutoff A are chosen to be 5 and 10000, respectively. The damping 
constant 7 is 0.5. 

some discretion is needed in the treatment of the cutoff scale. If we ignore this contribution, 
then one will encounter the following unphysical situation: If the oscillators are initially 
in a separable state , prepared as Gaussian wavepackets, their final state is always, at 
least marginally, entangled even though the mutual interaction is turned off. In contrast, 
if the cutoff contribution is accounted for, then the final state of the oscillators will not 
be entangled unless their mutual interaction is strong enough. Secondly, the cutoff scale 
always appears in the form In A, so the separability criterion is not very sensitive to the 
choice of the cutoff scale unless it takes some extreme values. In Fig 3, we let the cutoff 
scale A go up to a very high value relative to uj. We see that the separability criterion C _|_ 
become always positive above a critical value Ac, and Ac is highly sensitive to the choice 
of 7. Comparing the two plots in Fig 3, we see a mere change in 7 causes a dramatic shift 
in the value of Ac. Generally speaking, only for very weak oscillator-bath coupling will the 
cutoff-dependent terms play a subdominant role in the separability criterion. 

So far we have presented the general features in how the separability criterion depends 
on the interactions. We now investigate the role of temperature in the criterion. 

6.2.2 low temperature f3uj 3> 1 

Generally speaking, with increased temperature, thermal fluctuations will become increas¬ 
ingly important in affecting the dynamics of the oscillators from their respective baths. 
Quantum coherence is expected to deteriorate accordingly. We expect similar degradation 
may occur in entanglement. Thus it is reasonable to conjecture that once the temperatures 
of the baths are raised above a certain critical value, the degradation can be so severe that 
the oscillators become completely separable. However, the situation is more complicated 
for the present setup because two independent thermal baths are involved. It turns out 
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Figure 5. The symplectic eigenvalue ? 7 < is plotted with respect to the temperatures of two private 
baths. The black curve r]^ = 0 divides the separable state (ry< > 0, pink shade) from the entangled 
state ( 7 y< < 0, green shade). Essentially the curve traces along the region /3a; = 0{1). This result 
can be easily mapped to the logarithmic negativity by Ejs/{p) = maxjO, — ln 2 ? 7 <;}. The oscillating 
frequency oj and the cutoff A are chosen to be 5 and 10000, respectively. The damping constant 7 
is 0 . 1 , and the inter-oscillator coupling a is 20 . 


that lowering the temperature of one of the thermal baths does not necessarily guarantee 
entanglement between the oscillators. Thus the concept of a universal critical temperature 
is less well-defined in multiple bath situations. 

Here, we discuss the functional dependence of the separability criterion on tempera¬ 
tures. To begin with, let us suppose that it takes on a generic form C_|_ = C(/5i) h)- When a 
steady state is reached, the separability criterion should be invariant under the exchange of 
f3i and (32 because the configuration of the total system is designed to be symmetric when 
we swap one oscillator and its private bath with the other oscillator and its private bath. 
This implies that (+{(31,(32) = (+{(32, (3i). However, it is unlikely that the temperature 
dependence of the separability criterion can be reduced to a function of \(3i — (32\ solely. 
If C-i- were a function of \(3i — (32\, it would imply that the separability criterion could be 
independent of temperature for the case (3i = (32 where it would further suggest that both 
oscillators should be either separable or entangled for all temperatures. We have shown 
that at least they can not always be separable because in the zero temperature case, we 
found that both oscillators can be entangled for certain choices of parameters. On the 
other hand, it is hard to believe that both oscillators remain entangled even at very high 
temperature. Thus we rule out the possibility that the separability criterion may depend 
on |/3i — / 32 |- The same features are also shared by the symplectic eigenvalue rj+, as can be 
seen in Fig 5, but there are two differences: r/< is monotonic with respect to the parame¬ 
ters of the joint system and it does not rise up as steeply as the separability criterion in 
the high temperature regime. The latter is related to the extra factor — 1/4) in the 
criterion. Furthermore we observe that even for (3i 7 ^ (32 where the reduced system is de¬ 
scribed @by?@ asymmetric Gaussian states, the symplectic eigenvalue rj+, thus negativity, 
still gives a consistent and physical picture with respective to the ordering of the density 
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matrix in terms of the relevant parameters in question. 

Thus, to dehne more precisely a critical temperature /3c, we will look at the special 
case of /3i = / 32 . In this case both private reservoirs have the same temperature, yet they 
are totally uncorrelated. This setup is still distinct from the case that two oscillators share 
a common bath. In the shared bath case, the oscillators can influence each other indirectly 
through their interaction with the same bath, whereby non-Markovian effects enter in 
their dynamics, with dependence on their spatial separation (see, e.g., [9]) . This type 
of effects are absent in the private bath configuration; nonetheless, other than the direct 
influence from its own bath, each oscillator can still experience, by means of direct coupling 
between the two oscillators, the action of the other bath associated with the other oscillator. 
Therefore the equal-temperature private baths and the single common bath configurations 
are not the same, but, as we shall see, there are some similar features. Moreover, in this 
special case the two-mode Gaussian state becomes symmetric so the negativity can give an 
unambiguous comparison of entanglement between states. 

In the low temperature limit, we hnd the finite temperature correction to the separa¬ 
bility criterion is given by 
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where C,^'^ is the zero-temperature result in (6.13). It is interesting to note that the cor¬ 
rection may change sign as the inter-oscillator coupling a varies from zero to its upper 
limit. The upper limit of a is constrained by the condition Q- = = 0, so 

CTmax ~ O(w^). When a = 0, the term linear in 7 vanishes due to io+ = iO- there, but 
the term quadratic in 7 is positive. Hence the correction starts off with a positive value. 
On the other hand in the limit U- —)■ 7 (i.e., a —)• (Tmax), we hnd the hnite temperature 
correction gradually becomes negative 
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where we have used the L’Hopital’s rule to evaluate the limit of such an expression 

1 

1 In X T 

lim X In — = — lim ^j— = — lim .. = lim x = 0 . (6.19) 

X x^O -L J - 

X x 2 

This feature reveals the non-monotonicity of the separability criterion at hnite temperature. 
We stress that this errant behavior does not affect us from reading off the critical values 
of the parameters. 
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Figure 6 . The separability criterion is plotted against the inter-oscillator coupling cr and the 
oscillator-bath coupling 7 at low temperature. In each plot, we show the ('+ curve for three different 
bath temperatures. We see that the critical temperpareue is higher for stronger a but weaker 7 . 
The oscillating frequency lo and the cutoff A are chosen to be 5 and 10000, respectively 


Next we look into the effect of finite temperature correction on the critical value of ac 
where C-i- transits from a positive to a negative value. In the zero temperature case, we 
have found this critical value in (6.16), now denoted by Generally speaking the finite 
temperature correction of C-i- does not necessarily vanish at fi = as does. Instead 
we find at cr = cTc^^ the finite temperature correction of C-i- is 


27^ 

3/32u;4 


, A 

In- 1 

u> 


+ • • • , 


( 6 . 20 ) 


which is always positive. It means that this correction shifts the curve C+ upwards by 
about the order of magnitude 0 ( 7 ^). It thus implies that the critical values of a will 
increase because in general C-i- decreases with a, as seen in Fig 6 . In addition, a higher 
bath temperature results in a larger correction, and in turn causes ac to be even greater. 
Therefore thermal fluctuations from the baths make entanglement harder to maintain. The 
higher the bath temperature, the more severely the entanglement will deteriorate. This is 
consistent with our expectation. However, we may be concerned with a possible loophole 
related to what we found earlier that the finite temperature correction to C+ may change 
sign with increasing a. If it occurred before ac, we may encounter the opposite conclusion 
that the lower bath temperature will instead do more harm to the quantum entanglement 
in the system. We will argue that this is not the case. Let be the low-temperature 
correction, so that ()+ = C+^ • Since we have shown that when = 0, we have 

> 0, it implies that when = 0, we have C,+ '^ < 0 but remains positive. On the 
other hand, since monotonically decreases, = 0 will not occur until a > ac- Thus 
the separability criterion still offers consistent predictions on the behaviors of the critical 
parameters. 

From (6.17) we can derive a relation among the critical values of 7 , a and /3 for the 
small 7 cases. Similar to (6.15), we can show in the low temperature limit that 
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Figure 7. The trend of with respect to the inter-oscillator coupling a and the inverse temper¬ 
ature 13 when both private baths have the same temperature I3~^. It can be translated into the 
logarithmic negativity by Ej^{p) = maxjO, — ln277<}. We also draw a reference line 77 < = 1/2, the 
region below which represents the entangled final state of the joint system. In addition, all these 
curves are monotonic with respect to the parameters in discussion. The oscillating frequency w and 
the cutoff A are chosen to be 5 and 10000, respectively. The damping constant 7 is 0.5. 
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( 6 . 21 ) 


Note that the expression in the second pair of square brackets is negative in the low 
temperature case. We have argued earlier that at low temperature we don’t need a strong 
inter-oscillator coupling to safeguard quantum entanglement, so the curve C-i- can vanish 


for the small values of a. Furthermore if A 3> w, we have In 
safely conclude that 
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for small a and ui, and that the second pair of square brackets in (6.21) is negative in 
the low temperature case. Alternatively we may roughly see this based on the arguments 
that for the expansion to be valid we need cr/cj^ < 1 so the second term should be smaller 
in magnitude than the first term in the second pair of square brackets. The physical 
implication of (6.21) is that the critical temperature is lowered when the oscillator-bath 
interaction gets stronger and vice versa. 

Presently we have seen that the low-temperature correction can change sign for suf¬ 
ficiently large mutual coupling; however, this does not affect its usefulness to identify of 
the existence of entanglement. This unwelcome feature only makes murky the quantitative 
description of entanglement based on the separability criterion, and it can be traced back 
to the fact that the separability criterion contains not only but also rj^, whose existence 
distorts the information about entanglement, delivered by rj^. By comparison negativity 
is freed from this nonintuitive, unphysical behavior. Let us analyze the finite temperature 
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correction of t/<. In the same configuration, it takes a much simpler form 



15/34a;5^;5 


7 + 0(7^), 


( 6 . 22 ) 


with = (u^ihcr. We immediately see that it is always positive and monotonically increas¬ 
ing for all permissible values of a. Moreover, the finite temperature correction of ri<^ is a 
monotonically decreasing function of /3. This, unlike the separability criterion, give a plau¬ 
sible and physically intuitive description of the extent the state is entangled. Furthermore, 
since the analytical expression of r/< is much simpler than that of the separability criterion, 
it simplifies the analysis on the critical parameters. Expand out Uc = -!-••• 

and plug this expression back into the symplectic eigenvalue r/< = 77 ^^ -|- -!-••• = 1/2 

where satisfies = 1/2 at zero temperature. We find 
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which is positive-definite for all permissible ranges of the coupling constant a and tem¬ 
perature Note that the dependence on the cutoff scale is hidden in the expression of 

Finally we calculate the critical temperature via the criterion rj^ = 1/2 in the low 
temperature regime. Expanding with respect to large [3 gives 
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(6.24) 


where f{z) is defined in (4.5). Solving = 1/2 leads to 

P = _ ( 7 ^+ 7 <^)^ _ 


(6.25) 


The inverse critical temperature /3c grows with increasing 7 but falls off with increasing a. 

Therefore, we can see that in the low temperature regime the critical temperature jSc at 
which ? 7 < = 1/2 is higher for stronger inter-oscillator interaction, and for weaker oscillator- 
bath coupling 7 . This is totally in line with our intuition that the temperature and the 
oscillator-bath coupling 7 will corroborate to disrupt the quantum coherence between the 
oscillators and make them harder to remain entangled, while the inter-oscillator coupling 
will tend to enhance the coordination of both oscillators so their entanglement become 
more robust. 


6.2.3 high temperature /3a; <C 1 

We now turn our attention to the high temperature regime and ask if entanglement at high 
temperatures is at all possible. 

From the plot of the symplectic eigenvalue rj^ against the bath temperatures /3i and 
P 2 in Fig 5 we see that the surface rj^ forms a very flat basin which is symmetric with 
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Figure 8 . The symplectic 77 + is plotted against temperature. We show the numerical result 
and the low-, high-temperature approximations of ? 7 <. For stronger inter-oscillator coupling, the 
high-temperature approximation yields a very satisfactory result in the region where the transition 
occurs, in comparison with the numerical calculations. The oscillating frequency w and the cutoff 
A are chosen to be 5 and 10000, respectively. The damping constant 7 and the inter-oscillator 
coupling a are 0.2 and 24 respectively. 


respect to /3i and / 32 . The surface r]^ will mildly rise up in the vicinity of Pi ^ 2 ^ = 0(1) 
when we approach from the low temperature end. Next we sec it sharply climbs up, 
crossing the dividing curve = 1/2 in the region = 0 ( 1 ), and enters the high 

temperature regime. Thus we can make a first observation that, roughly speaking, the 
dividing curve of r]^ = 1/2 follows Piuj = 0(1) and then turns to ^ 2 ^ = Secondly 

it implies that separability is determined by the temperature of the warmer bath, instead 
of the temperature difference, as was mentioned in the previous section. Thirdly, since 
from earlier discussion we know entanglement tends to survive at higher temperature if the 
mutual coupling between oscillators is stronger, we use the high temperature approximation 
to find the critical temperature in the strong a regime. As shown in Fig 8 , we compare 
the numerical calculation of ?]< with its low and high temperature approximations, and 
see that in the large a case the high-temperature approximation yields a very consistent 
behavior of rj^ in the vicinity of /3a; ~ 0{1), in comparison with the numerical results. 

In the high temperature approximation, the symplectic eigenvalue rj^ is given by 
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(6.26) 


The cutoff-dependent factor in those higher order expressions is less important in the weak 
7 limit because it always appears with the small parameter 7 / 0 ;. The critical temperature 
occurs at Directly solving a transcendental equation like (6.26) for /3c is not 

possible. Nonetheless since In /3A always pairs up with 7 , we can use the iteration method 
to derive /3c. If we first ignore terms of the order 0 ( 7 ) and higher, we find /3c is given by 
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Figure 9. We plot the separability criterion C_|_ and the symplectic eigenvalue ? 7 < together. They 
crisscross at the critical temperature, and therefor give the same information about the existence of 
entanglement. However, separability criterion falls off and rises up with increasing 13. Note that we 
shift the values of downward by 1/2, That is, what we plot in fact is ? 7 < — 1/2. The oscillating 
frequency w, the cutoff A, the damping constant 7 and the inter-oscillator coupling cr are chosen to 
be 5, 10000, 0.2 and 24 respectively. 


2\/3/\/3a;2 -|- 4a. Substituting it back to seek a correction of order £>( 7 ) we obtain 
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It is indeed consistent with the statement that jScUJ = 0(1), and it rules out the possibility 
of the existence of entanglement in the regime j3u <C 1. Again it reveals the fact that with 
larger inter-oscillator coupling a we see a higher critical temperature; on the other hand, 
stronger oscillator-bath interaction 7 will cause the critical temperature to decrease. 

The same results can be found if we investigate the high temperature approximation 
of the separability criterion. This is no surprise since we have previously discussed that 
the separability criterion is perfectly valid for identification of entanglement except for a 
quantitative measure of entanglement. For example, as shown in Fig. 9, the separability 
criterion C-j. and the symplectic eigenvalue 7 < give the same prediction about the location 
of the critical temperature, but the separability criterion is not a monotonic function of 
the temperature, which makes it inappropriate as an entanglement measure. 

With temperature measured in ratio to the oscillator’s natural frequency /3cu we can 
conclude that there is no high temperature entanglement in Case Cl, namely, between two 
oscillators each interacting with its own bath. 
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7 Intuitive Understanding of Entanglement Behavior 


So far we have taken quite some labor to show that asymptotic entanglement between 
oscillators are easier to sustain for stronger inter-oscillator coupling but weaker oscillator- 
bath interaction. 

Here we would like to offer a physically more transparent illustration as to the compet¬ 
ing roles between these two kinds of interactions in terms of normal modes of the oscillator. 
The Langevin equations (2.1) and (2.2) can be easily decoupled by forming two normal 
modes 


X+ = 


Xl +X2 


X- = Xi - X2 , 


(7.1) 


and the corresponding equations of motion are 


X+ + ‘^lX++^Xx+ = ^(^ 1 +^ 2 ) , 

X- + 27x_ -k X- = — (u - 6) • 

m ^ 


(7.2) 

(7.3) 


Since we are interested in the late-time dynamics, we will not write down the homogeneous 
solutions to the Langevin equations. Following the earlier discussions we find that the 
inhomogeneous solutions are given by 


X-h ^ + ^2 (s) 
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(7.4) 

(7.5) 


with 


^ '^’^sinH±r, — 7 ^ . (7.6) 

The frequency is the resonance frequency of the normal modes x± ■ Hence the stronger 
inter-oscillator coupling a implies smaller values of a;_ but larger values of a;+, and in 
turn smaller and larger Since the amplitude of the normal modes x± is related to 
the ratio 7 /H±, stronger inter-oscillator interaction will induce a larger amplitude of the 
mode x_, which will grow with increasing values of cr, meanwhile it causes the mode x+ 
to oscillate subdominantly and its amplitude decreases with a. This is intuitively plausible 
since, e.g., for a very soft spring, or for a particle in a very shallow harmonic potential, 
a small disturbance could easily induce a large displacement in its motion. Thus in these 
circumstances it tends to have a large position uncertainty. Furthermore, the consequence 
from the normal-mode dynamics hints at the fact that when we form the displacements of 
two oscillators by superposing the normal modes 

Xl = X-H + ^ X- , X 2 = X-H - ^ X- , (7.7) 

the mode X-i- can be overshadowed by X-- The original displacements xi; X 2 are more 
or less determined solely by the mode X-j ia particular in the strong mutual coupling 
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limit —)• 0. Furthermore, in this limit, xi X 2 will be out of phase by vr. Likewise, 

following similar arguments and taking care of contributions from the resonance, we see 
that in the case of the conjugate momentum pi, p 2 of the two oscillators, the contribution 
of can dominate over that of p- for strong mutual coupling between the oscillators. 
Furthermore, the phase difference between xi X 2 is reflected by the fact that in this 
limit we should have Vn ~ — V13. It is particularly easy to see this for the special case 
/?! = (32 ■ The formal late-time expressions of Vn and V13 in this case are 



in the limit a ^ uP' where |c^ > \^ 2 ^\f<-)\‘^■ In addition Eqs. (5.1) and (5.3) also ex¬ 

plicitly demonstrate the same relation. Similarly p_|_ is the dominant mode in the conjugate 
momenta pi, p 2 , so we may expect V22 ~ + V24 and this is clear from 



for the case /3i = /32. 

Now let us take a look at the formal expression of the symplectic eigenvalue 7/<. From 
(6.10), we see 

V< = ({xi,X-h})({vi,V-}) = |({x+,X-h})({v-,V-}), (7.8) 

due to the fact that there is no cross-correlation between two normal modes in the case 
/3i = (32- It is clearly seen that r^< is composed of subdominant contributions only, which 
all have smaller uncertainty with larger a. Moreover they decrease with increasing values 
of the inter-oscillator coupling a. This makes possible that the symplectic eigenvalue r/< 
can fall off with strong inter-oscillator coupling, thus that entanglement can be sustained 
at higher temperature. 

In summary we show that by analyzing the behaviors of the normal mode frequencies 
with respect to various couplings and parameters of interest, we may get an intuitive 
understanding of the general features of the entanglement between two oscillators in relation 
to these parameters. 


8 Summary of Results and Comparisons 

8.1 Summary of Results for Entanglement in Systems in NESS 

Having shown the quantitative details in the last section we now provide a summary of 
the qualitative features of entanglement dynamics in the case (Case Cl) studied here for 
quantum systems in NESS. For two bilinearly coupled oscillators each interacting with its 
own bath, we find: 
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1. Quantum entanglement in systems of this setup is harder to sustain at finite tem¬ 
peratures. Thermal fluctuations from the baths disrrupt the coherence between the 
oscillators. 

2. Both the separability criterion and the negativity are perfectly good indicators for 
the existence of entanglement. However, the former is not necessarily a monotonic 
function of the parameters in the configuration, so it does not qualify as an entangle¬ 
ment measure. It cannot give a consistent, quantitative comparison between different 
entangled configurations. 

3. The entanglement criterion or the symplectic eigenvalue ??< in general is not a 
function of the bath temperature difference, in contrast to thermal transport in the 
same setting [ 1 ]. 

• Lowering the temperature of one of the thermal baths does not necessarily help 
to keep the entanglement between the oscillators. 

• The notion of a critical temperature, where C-i- = 0 or r 7 < = 1/2, is better defined 
when two private baths have the same temperature. 

4. Entanglement between the two oscillators is reduced for stronger oscillator-bath in¬ 
teraction, but enhanced for larger inter-oscillator coupling. They play competing 
roles as far as their influence on entanglement is concerned. 

• strong inter-oscillator coupling better links the dynamics of the two oscillators, 
and thus improves the coherence between them. 

• uncorrelated environment fluctuations corrupts the correlations between the os¬ 
cillators; stronger oscillator-bath interaction will compound this effect. 

5. For weak oscillator-bath coupling the critical temperature satisfies ficOJ rsj 2(1 + 
4cj/3a;^) ^. This supports the rough estimate condition (3cUJ 0 ( 1 ). 

• For strong oscillator-bath coupling the critical temperature depends on the 
damping constant 7 and the environment cutoff frequency A. 

• The effect of environment cutoff cannot be ignored in the low temperature and 
the strong oscillator-bath coupling cases. 

6 . Asymptotic quantum entanglement disappears in the high temperature regime /3a; <C 
1. There is no hot entanglement in systems (with bilinear constant coupling) under 
NESS conditions. 

8.2 Comparison: System in a Private Bath vs in a Common Bath 

It is useful to make a comparison of the case studied here (Case Cl) with what we have 
studied in Paper I (Case B), namely, a system of two bilinearly coupled oscillators inter¬ 
acting with one common bath. 
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Case B: common bath. In Paper I we have studied the case of two coupled oscillators 
at a finite spatial separation, both interacting with a common thermal field bath, which 
is a finite temperature generalization of the work [9, 30]. For comparison with Case Cl 
studied here we only need to consider the limiting case when the two oscillators are placed 
next to each other thus ignoring the spatial variation of entanglement. The action for 
this setup takes the form 
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Since the two oscillators share the same bath, we can decompose the two harmonic oscillator 
variables into the fast mode (or center of mass) and the slow mode (or relative coordinate) 
variables, x+ = \{xi + X 2 ), X- = Xi ~ X 2 whence the action becomes 


S = 


dt\ 


m 


muj. 


X- 


- ,,2 


T- 


+ 2 dt 


m 


X+ 


mu). 


+ ...2 


x+ 


+ 2e i drx x+{t)^{x — z{t)) (l){x) + d‘^x - d4>{x) , (8.2) 


1 2 


where ± u. 

We see that the slow mode x- is decoupled from the bath, while the fast mode x+ now 
interacts with an effective bath, described by the same scalar field but with the reduced 
amplitude, cj)jy/2 and with an effective coupling strength enhanced to \/2e. The Langevin 
equations for the fast and slow variables are 

X+ + 27 X'+ + , (8.3) 

m 

X-+a;^ = 0. (8.4) 


Note there is a subtle issue in this pair of Langevin equations. Although we still write 
the oscillating frequencies of the fast and the slow mode as u}±, now they have quite 
different physical contents. We observe that the fast mode couples with the bath, so its 
oscillating frequency will acquire a correction due to its interaction with the bath. This 
correction is absent for the slow mode. Nonetheless because the correction depends on the 
environment cutoff and it is of the order buP' ~ ©(yA/w^), we expect that in the weak 
oscillator-bath coupling limit, this correction is moderate. In fact it has been shown in [9] 
that the oscillator-bath coupling constant should be reasonably small or else in may induce 
instability due to the Coulomb-like interaction at short distances. 

The stochastic force ^ in this case still possesses the statistical properties 

(C(t)) = 0, = (8.5) 

^See discussions in [9] in how dose the two oscillators can be placed to avoid possible singular retardation 
effect and non-Markovian behavior. 
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The solutions to (8.3) and (8.4) are 


X+{t) = X+(0) + x+(0) + — [ ds d^ 2 ^\t - s)^(s), (8.6) 

m Jo 

X-(t) = d[~\t) x-(0) + X-(0), (8.7) 


with 


0 i+)(t) = e-T'* 


cos W+t + sin W+t 


0^ \t) = COS W-t , 


05 +^(t) = :^e-^*sinTT+t, Wl=ul-x‘^. 

( 8 . 8 ) 

t)(-) {t) = ^ smW-t, = wi . 

W- 

(8.9) 


Here we note that the slow mode in the common bath case is non-decaying, so some of the 
initial information of the system can be kept to the very end of the evolution. This is in 
strong contrast with the private bath case studied here. The corresponding component of 
the fast mode decays with time. Thus at late times the fast mode only responds to the 
environment. 

The original oscillator variables thus evolve according to 




2 L ^ 


'm Jo 
m Jo 


ds - s)Cis ), 

(8.10) 

ds - s)^{s), 

(8.11) 


at late time t S> 7 We find the elements of the covariance matrix in this case are given 
by 


ixlit)) 

ixlit)) 

^({Xl(t),X 2 (t)}) 

iplit)) 

ipUt)) 
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ixlit)) 
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9‘r’7i){xii(o)>+».r'"(‘)(x'-(o)) 

2 nt 




+ "W / dsds'di^\s)d^ 2 ^\s')GH{s - s'), 


m? Jo 
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(x“(0)) + 8r'"(f) (xii(0)> 


^(-) 2 m /V .2 


+ I ds ds'd^\s)i)^^\s')GHis — s') 


iplit)), 
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2 r 


-Wl 


(x“(0)) + 8r'"(f) (xii(0)> 


^.(-) 2 m /V .2 


( 8 . 12 ) 

(8.13) 

(8.14) 


(8.15) 

(8.16) 
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5 ({Xl((),P 2 (*))) 

\{{x2{t),nm 

^{{X 2 (*).Pl(*))) 


+ e^ 


ds ds' — s'). 


, m 

+^4 

m 

1 


v[-\t)i>[-\t) ixim) 
ixim 


l{{Xl{t),P2{t)}) . 


ixim 

+di-^\t)di-\t) ixim 


(8.17) 

(8.18) 

(8.19) 

( 8 . 20 ) 
( 8 . 21 ) 


For the last four elements, the term caused by the environment vanishes in the limit t —)• oo. 
The parameter A is used as a marker for the intrinsic components, and can be set to unity 
with no consequence. This is in contrast to the components induced by the environment 
which have dependence. 

One feature in the common bath case noticeably different from the private bath case 
is that the elements of the covariance matrix still contain the information about the initial 
conditions even though the system has evolved to late time. This is a consequence of the 
fact that one of the normal modes is completely decoupled from the bath such that part of 
the initial information is retained in the system and is not lost into the environment. On 
the contrary, for the private bath case, both the slow and the fast modes are coupled to the 
environment, and it causes the dispersion of the initial information into the environment. 
We also note that the components induced by the environment are typically smaller by an 
order 0{'y) than the components intrinsic to the quantum evolution of the oscillators. 

Moreover, we have shown that in the private bath case, stronger inter-oscillator cou¬ 
pling renders the oscillating frequency of the slow mode smaller than that of the fast mode. 
It implies that when the interaction between the oscillators are comparable with the orig¬ 
inal oscillating frequency co the slow mode will dominate over the fast mode. From (7.7), 
we see that the late-time correlation between xi and X 2 is prone to be negative, meaning 
that xi tends to be anti-correlated with X 2 - This is not the case for the common bath 
case. When both oscillators share a common bath, only the fast mode couples with the 
bath. The coupling between the oscillators plays a minor role because is always of the 
order cj. At late time t 7 “^, we see that both xi and X 2 are more or less led by the 
fast mode, apart from the intrinsic quantum evolution of the system inherited in the slow 
mode. Hence the bath tends to drive two neighboring oscillators into correlation, meaning 
that the correlation between xi and X 2 induced by the shared bath tends to be positive. 

At this point, some discretion is advised. First, we observe from (8.14) that the 
correlation caused by the intrinsic quantum evolution of the system carries a negative sign, 
and will counteract with the correlation induced by the environment, so the total correlation 
between xi and X 2 may not always be positive definite at late time. This will also introduce 
additional complexity in the entanglement for the shared bath case. Secondly, unlike the 
private bath case where the elements of the covariance matrix approach a time-independent 
constant when the NESS is reached, the corresponding elements in the shared bath case 
remain oscillatory in time reflective of the intrinsic quantum dynamics of the system. 
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We end with a summary of the qualitative behavior of entanglement in a system of 
two coupled oscillators comparing between the two cases; in one case this system interacts 
with a common bath, and in the other case, with their own private baths, but kept at the 
same temperature. 

1. From the structure of the normal modes, we find that 

• private bath; both degrees of freedom are coupled to the bath, so they behave 
like a pair of damped driven oscillators with different oscillating frequency. 

• common bath; only the fast mode is coupled to the bath, the slow mode is 
totally decoupled from the bath and it acts like a free oscillator. 

2. If we separate the elements of the covariance matrix into the intrinsic and the induced 
components; the former depends on the initial conditions of the oscillators but the 
latter is entirely driven by the environment, independent of the initial conditions of 
the oscillators. We see that 

• for the mode coupled to the bath, the intrinsic component will be exponentially 
small at late time, and only the induced component survives. 

• for the free mode, the intrinsic component oscillates all the time and there is no 
induced component. 

3. This implies that the initial Gaussian conditions will be washed out for the mode 
coupled to the bath, but they will survive at late times for the uncoupled mode. 

• private bath; the initial Gaussian conditions will be irrelevant to the asymptotic 
entanglement, 

• common bath; they remain significant, so the final state of entanglement de¬ 
pends on the choice of the initial conditions. 

4. At late times the entanglement measure for the private bath case is time-independent, 
but for the common bath it continues oscillating in time. 

5. The amplitude of the driven mode is related to the mode frequency. The smaller the 
frequency is, the larger the driven amplitude will be. 

• private bath; slow mode will have larger driven amplitude than that of the 
fast mode, so the dynamics of the original canonical variables, which are the 
superposition of these two modes, will be dominated by the slow mode, especially 
when the mutual interaction is strong. 

• common bath; since there is one driven mode and it is the fast mode, the driven 
amplitude does not change too much as the mutual coupling strength varies. 
However, the asymptotic dynamics is determined by the relative magnitude 
between the slow mode (intrinsic component only) and the fast mode (induced 
component only). 
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— if the fast mode dominates, then the asymptotic elements of the covariance 
matrix will be more or less constant in time with small ripples. 

— if the slow mode dominates, then they will oscillate in time. 

6. The inter-oscillator coupling (cr > 0) plays a more important role in the private bath 
case, but a minor role in the shared bath case. 

7. In the private bath case, entanglement is easier to survive for stronger inter-oscillator 
and weak oscillator-bath coupling, but in the shared bath case, both factors can be 
overshadowed by the intrinsic components, which are sensitive to the initial conditions 
of the oscillators. 

8. The asymptotic entanglement criterion in the common bath case can thus be broken 
into three components: one involving the fast mode only, one with slow mode only 
and the cross term. 

• if the fast-mode part is subdominant, then the resulting entanglement criterion 
will oscillate with time, and that can cause sudden death [28] and revival [29-31] 
(SDR). 

• if the fast-mode part is dominant, then there is no SDR phase. 
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A Expressions for Vis^t), Vi 2 {t), Vu{t), 1^2(t), 1 ^ 32 (t), V24(t) 

A.l Visit) = {{xi{t),X2{t)})/2: 

From (2.18), we see that 


Vi 3 = lim FL 3 (t) 

r^oo 



dK 
— nn 27T 


* (K)Df (k) («) + * (k)D 


)f( 


G 22 
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(k) 


dK 0 -(k 2 -a; 2 )[G]yL(K) + Gf (k)] 

y_^27r [(K2_a;2)2+4^2^2][(^2_^2)2+4^2^2] ’ ^ 


where 


D“*(k)D21(k) 

I)f*iK)&^iK) 


a [(k^ — cj^) — i 27^] 

[{k? - a;2 )2 + 472^2] [{k? - a;2 )2 + 472^2] ’ 
a [(k^ — + i 2 xk\ 

[(k 2 — a;^)2 + 472^2] |^(k 2 _ ^2 )2 _j_ 4-j/2^2j 


We see that both baths contribute equally. 


(A. 2 ) 

(A.3) 


A.2 Vi 2 it) = {{xi{t),piit)])/2-. 

We find 

\ ({xiW. Kit)}} = mD;‘(t)Df (i) {xlm + 1 (pl(0)) 

^ III 

p 2 /■* 

+ — / dsds' - s)Df (t - s') G^^(s - s'). (A.4) 

4o 


Note that the overdot represents time derivative with respect to the argument of the 
variable. The late time limit of Vi 2 (i) is then given by 


V 12 = lim Fi 2 (t) 

t^OO 


^ d^ 

m y_oo 27r 


Knl'^*iK)-Dl>^iK) G’k^iK) = 0. 


(A.5) 


The result is identically zero because the integrand is odd in k. 


A.3 Vu{t) = {{xiit),P2it)})/2: 

The late-time limit of Vuit) can be inferred from (A.4), 

2 r°° Pk ~ ~ 

Vll = -i- —KJ^f*iK)Gf{K)G^j^iK) 
m 2 -k 

_ 1'°° dK 2axK‘^[G^^{K) - Gf^iK)] 

m 27r [(k^ — CU^)^ + 472 K 2 j |^(^2 _ ^2 ^2 _j_ 4 ^ 2 ^ 2 j 


It is interesting to compare this result with (A.l). The latter depends on additive contri¬ 
bution from both baths while the former has a subtractive contribution between baths. 










A.4 V22{t) = {{pi{t),pi{t)})/2 = {pl{t)): 

By a similar derivation of ({xi(t),Xj(^)})/ 2 , we find 

^{{Pi{t),Pj{t)}) = (t)Df (t) ixliO)) + D*2^(t)D^^(t) {pl{0)) 

+ e^ [ ds ds' t)f{t - s)Df (t - s') G^^(s - s'). (A. 7 ) 

Jo 

Thus we have the late-time limit of V22(t) given by 


/ OO 

-oo 



This is similar to Vn. 


ds ds' D^^(s)Df (s') Gfis - s') 
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|2 r^22/ 


(A. 8 ) 


A.5 V32{t) = {{ x 2 { t ), Pi { t )})/ 2 : 

We show that V32(t) = — Ti4(t), that is 


V32 = —i 


dK 


dn 


= +- ir- 


D2i*(Ac)DinAc) G^^{k) + G“(k) 

2axK^[G]^{K) - Gf (k)] 
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= -Vi4. 


(A. 9 ) 


In the context of nonequilibrium transport, Vi4(^) is related to the power done by the 
oscillator 1 to oscillator 2 by means of the mutual coupling, while V32(0 is the other way 
around. At late time when the steady energy current is established, both should be equal 
in magnitude but opposite in sign. See. e.g., [ 1 ] 


A.6 V'24(t) = ({pi(t),P2(t)})/2: 

Eq. (A. 7 ) tells us that at late time t —>■ 00, the element V2i{t) becomes 
V24 = e^ G“(^) 

J — OO 

_ 2 f°° dK crK^(«;^ - w^) [G}^(k)-I-G|^(k)] 

27 r [(^2 _ u ;^)2 -I- 472^2] [(k 2 — u ;^)2 + 472^2] 

Again it has a similar structure to V13 in (A.l). 


B High Temperature forms of V 13 , V 14 , V 225 V 24 

At high temperature the Hadamard function G^(k) is approximately given by 

= (B.l) 
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From (A.l) we need the integral 
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47 — ( t 2 ’ 

to evaluate V13 at high temperature. Thus we have the high temperature limit of V13 given 
by 
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27 
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1 1 
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(B.3) 


2m oj^ — 0-2 Y/3i (32, 

It means that in this configuration, xi X 2 anti-correlated and this anti-correlation 
grows with the mutual coupling strength a. 

To calculate V14 we need the integral 

..2 ^ 1 
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TT 


[( k 2 - ^2)2 + 4^2^2] [(^2 _ ^2 )2 + 4^2^2] 4^ 4^2^2 + 


cr^ 


so that from (A. 6 ) in the high temperature limit, V 14 is given by 

4720- [1 ll f 1 1 

Vl4 =- 24 


TT 
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4W272 4 - (j 2 \ 


(B.4) 


(B.5) 


The correlation between xi and p 2 diminishes with increasing mutual coupling. Moreover, 
the correlation disappear when both thermal baths have the same temperature. 

Similar to those in evaluating Vn, the following two integrals are needed for evaluation 
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The high temperature form of V 22 = (Pi(c>o)) is given by 
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In the case /3i = /S = / 32 , h reduces to 




(B. 8 ) 


(B.9) 


which is independent of both coupling strengths 7 and a. 
Finally for V24 the integral 
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(B.IO) 


vanishes identically, so V 24 = 0 in the high temperature limit. Thus we need the next order 
contribution. If we expand the Hadamard function in (B.l) one more order in /3, we find 
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(B.ll) 
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Then we need the integral 
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and from (A. 10 ) we obtain 
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This contribution is relatively small in the high temperature limit /?a; —>■ 0. 


(B.13) 


C Zero-temperature expressions for V13, V14, V22? V24 
C.l Vi3 

We first evaluate the integral 
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We thus obtain V13 


Vg' = ^J3 = — 
vrm 8 m 


f{n+) fin.) 


n. 
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C .2 Vi 4 


The elements vanishes because the contributions from both thermal baths cancel. 

C.3 V 22 

Here comes the tricky part because divergence emerges when we evaluate the vacuum 
component of V 22 - We first calculate the following two integrals 
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Apparently J5 diverges since its integrand behaves like The sum of J5 + Jg is given 
by 
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The logarithmic divergence is regularized by a frequency cutoff A. Therefore the vacuum 
contribution of V 22 is 
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Now let us check some limiting cases of (C. 6 ), 
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1 . IT — 7 - 0 : when the mutual coupling is vanishingly small, the momentum uncertainty 
of Oscillator 1 becomes 

=-— In-+ —- ^ - f{n), (C.7) 

IT —>-0 TT OJ 4 sZ 

where 0^ = 7 ^ is the resonance frequency of the oscillator. This is the momentum 

uncertainty of the uncoupled oscillator when it couples to the vacuum fluctuations of 
the bath. 

2 . 7 — 7 - 0 : the leading contribution of the momentum uncertainty at late time in this 
case is 

"^22 = J (^+ + = ^ + C’(7) • (C. 8 ) 


C.4 V 24 
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Thus we have the vacuum component of V 24 given by 
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D Low temperature correction expressions for V13, V14, V22? V2 
D.l Vi3 


We first evaluate the integral 
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This implies that the low temperature correction to V 13 is given by 
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For Vi 4 we need the integral 
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Here the finite temperature correction behaves like so we will acquire a factor 
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once we consider all algebraically equivalent contributions in (4.7). Thus the correction to 
Vi 4 is given by 
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D.3 V 22 

Before evaluating V 22 ) we first evaluate the following two integrals 
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We have the low temperature correction to V 22 given by 
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The corresponding finite temperature correction to the kinetic energy of Oscillator 1 is 
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Here we need the integral 
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Therefore V 24 becomes 
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